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ABSTRACT 

A Bayesian inversion method is introduced for retrieving the fraction of ground flashes in a set of N 
lightning observed by a satellite lightning imager (such as the Geostationary Lightning Mapper, GLM). 
An “exponential model” is applied as a physically reasonable constraint to describe the measured 
lightning optical parameter distributions. Population statistics (i.e., the mean and variance) are invoked to 
add additional constraints to the retrieval process. The Maximum A Posteriori (MAP) solution is 
employed. The approach is tested by performing simulated retrievals, and retrieval error statistics are 
provided. The approach is feasible for N> 2000, and retrieval errors decrease as N is increased. 
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1. INTRODUCTION 

A constrained mixed exponential distribution model and a Bayesian inversion method are introduced for 
retrieving the fraction of ground flashes in a set of flashes observed by the future Geostationary Lightning 
Mapper (GLM); the method can also be applied to low Earth orbiting lightning imagers. The approach is tested 
by performing simulated retrievals, and retrieval error statistics are provided. 


2. DISTRIBUTIONS OF A LIGHTNING OPTICAL PARAMETER (SHIFTED MGA) 

Figure 1 shows the ground and cloud flash distributions of y = x - 64 km 2 (where x = the Maximum Group Area, 
MGA, in a flash). These data were derived from an analysis of 5 years of OTD data (Koshak, 2010), and 
discussed in Koshak and Solakiewicz (2010). The red analytic curve in each plot is an exponential distribution 
with a mean equivalent to the data average shown in the upper right comer of the plot. 
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Figure 1. Distributions of the (shifted) MGA variable, y. 


3. RETRIEVAL METHOD 

For a mixture of ground and cloud flashes, one can consider a superposition of exponential distributions. This 
gives, in general, the following Mixed Exponential Distribution Model, and Bayesian Inversion scheme: 


Mixed Exponential Distribution Model 

Basic Definitions: 

x - Maximum Group Area (MGA) in a flash 
y = x-64km 2 (Shifted MGA) 

N g = # Ground flashes in lat/lon bin of interest 

N-U Flashes in lat/lon bin of interest 
a - N z / N (Ground flash fraction) 

// = Population mean y for ground flashes (in lat/lon bin of interest) 
pL c - Population mean y for cloud flashes (in lat/lon bin of interest) 


Distribution of MGA modeled as a Mixed Exponential Distribution: 
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Bayesian Inversion 


Bayes’ Law: 
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Find parameters v- (a,p g ,p c ) that maximize the probability on LHS. 

This means one maximizes the following : 

m m 

S(v) = ln[ P(y |v)^(v)] = lnn J P(.v,|v) + ln *(▼) =Z ln 

M M 

Formally : 

= 0 v = "Maximum A Posteriori (MAP) Solution" 
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Practically : 

Use Broyden-Fletcher-Goldfarb-Shannon variant of Davidon-Fletcher-Powell numerical method to 
minimize - S(\). Also, R(v) is simplified by assuming model parameter independence, with P(a) 

uniform, and R(/z ) & R(/z ) both normal distributions. 

£ £ 


« -y,i» r , ( 1 -<*)_- 


+ 


Jit A 


M 


g 


Me 


+ In P(v) , 


4. GROUND FLASH FRACTION RETRIEVAL ERRORS FROM SIMULATIONS 


To test the retrieval method, lightning sources were simulated; the known ground flash fraction was 
compared to the retrieved value. A few examples of retrieval errors are shown in Figure 2 for TV=2000. 
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KEY 

[Ground Flash Fraction (alpha) Retrieval Error] 

Low Error: Green Shades 
Mid Error: Blue Shades 

High Error: Red Shade 



Note: All errors decrease when 
N increases. 
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Figure 2. Mean ground flash fraction (alpha) retrieval errors as a function of the known population 
mean values of MGA (2 left plots); the true alpha range is given at top of each plot. The standard 
deviations about the mean are shown in parentheses within each cell. Variation of retrieval error is also 
shown as a function of entire true alpha range (right-most plot). 


5. GLOBAL RETRIEVAL OF GROUND FLASH FRACTION (FROM OTD DATA; Figure 3) 

The ability to retrieve ground flash fraction has important benefits to the atmospheric chemistry community. For 
example, using the method to partition the existing OTD/LIS satellite global lightning climatology into separate 
ground and cloud flash climatologies would improve estimates of global lightning nitrogen oxides (NOx) 
production; this in turn would improve both regional air quality and global chemistry/climate model predictions. 
A first attempt is made here to retrieve the ground flash fraction on a global scale using just OTD data. 



For the 562 4x4 degree lat/lon bins shown 
to the left (each bin has N >2000), the 
following stats were obtained: 


Mean Alpha =0.151 
Std Dev. Of Alpha = 0.081 
Max Alpha = 0.739 
Min Alpha = 0.019 
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Figure 3. Ground flash fraction on a global scale using Optical Transient Detector (OTD) data. 
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